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We discuss the Schrodinger functional in lattice QCD with staggered fermions and relate it, in the classical 
continuum limit, to the Schrodinger functional regularized with Wilson fermions. We compute the strong coupling 
constant defined via the Schrodinger functional with staggered fermions at one loop and show that it agrees with 
the continuum running coupling constant in the Schrodinger functional formalism. We compute this running 
coupling in the "weak coupling phase" of many flavor QCD numerically at several values of the bare coupling 
and for several system sizes from L/a = 4 to 12. The results indicate that the /3-function for 16 flavors has the 
opposite sign than for few flavor QCD, in agreement with a recent claim, and with the perturbative prediction. 



1. THE SCHRODINGER FUNCTIONAL 
WITH STAGGERED FERMIONS 

The Schrodinger functional describes the evo- 
lution of a state at (Euclidean) time t = to 
another state at time t = T . Using the transfer 
matrix it can be written as a path integral with 
fixed boundary conditions at time t — and T . 
For staggered fermions all degrees of freedom can 
be fixed at both boundaries [|J. The Schrodinger 
functional can thus be represented as the path 
integral 



Z[W,(,t;W',('X'] = / [DU] 



(1) 



T-l 



n n Mx(^ 4 )cw^4)]e- 5G - ssF . 

Here W and W represent the boundary values 
of the gauge fields and £, Q, Q' and (' those of 
the fermion fields. The pure gauge action and 
measure are as in (2) . The fermionic part of the 
action is the usual staggered action with x an d X 
set to zero for x^ < and a; 4 > T, and to their 
boundary values at the boundaries |l],|| . 

Note that for staggered fermions the total num- 
ber of time-slices has to be even. Therefore 
the time extent T has to be odd. In the spa- 
tial direction we take the lattice to be of size 
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L (even!) and impose the generalized periodic 
boundary conditions [|| x( x + Lk) = e l9k x{x), 
X(x + Lk) =x(x)e- i6k . 

It can be shown (3) that for massless stag- 
gered quarks the classical continuum limit of the 
Schrodinger functional agrees with that from Wil- 
son fermions j^]. 



THE SCHRODINGER 
AT ONE LOOP 



COUPLING 



The "Schrodinger functional coupling con- 
stant" is defined from the response to a constant 
chromoelectric background field introduced via 
the boundary fields. We use those given in 
which depend on a parameter, 77, and define 
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with T 
scale, g 2 



L such that g 2 depends only on one 
g 2 {L). The normalization k has been 
chosen such that g equals the bare coupling at 
tree-level without any cutoff effects. 

The one-loop contribution from the staggered 
fermions to the coupling constant, g 2 = g 2 , + 
(pi,o + n /Pi,i)ffo + C(5o)j comes from the deriva- 
tive of the fermion fluctuation determinant 



Pis = 
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Table 1 

The first two 'non-log' terms in the expansion 
eq. (|]) of rj ^ is the average of the choices 
T = L±a. 
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-0.004416(1) 


0.00947(1) 




tt/5 


-0.00579695(2) 


0.009477(5) 




1.0 


-0.0068642(1) 


0.00947(1) 





Here rif = 4 for one flavor of staggered fermions, 
since they correspond to four flavors of continuum 
fermions. The fermion boundary fields are set to 
zero, and M is the fermion matrix. 

To be precise, the coupling is defined in [^Jlj as 
eq. (|J) with T — L so that it depends on a single 
scale, L. Unfortunately this is not possible for 
staggered fermions, since, as we have seen, L/a 
must be even but T/a must be odd. Instead, we 
average the couplings obtained with T = L + a 
and T = L — a to avoid an additional 0(a) effect. 

We have evaluated || of eq. @ for L/a 
ranging from 4 to 64, in steps of 2. Here the 
superscript ± stands for the choices T — L ± a. 
One can then extract the first few coefficients of 
the expected asymptotic form 

Pi,i{L/a)=r Q + s Q \og(L/a) (4) 
+ (n + Sl log(L/a))(a/L) + ■ ■ ■ . 

s should just be 2&o,i = — l/(127r 2 ), the 
fermionic contribution, per flavor, to the f3- 
function Q, and thus absorbed by renormaliza- 
tion. We indeed found this result. s\ was found 
to be compatible with zero. ro,i are listed in Ta- 
bic [l]. ro obtained from the choices T — L ± a 
always agreed to the accuracy given. 

r of eq. (Q) is the finite part of the fermionic 
contribution to the one-loop relation between 
bare lattice and running Schrodinger functional 
coupling. When converted to a relation between 
the MS and the Schrodinger functional coupling, 
we found || good agreement with the results 
of Sint and Sommer from Wilson fermions ||]. 
This confirms that our Schrodinger functional for 
massless staggered fermions is a correct regular- 
ization of the continuum Schrodinger functional 
with massless fermions. 
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Figure 1. The ratio 5 hl (a/L) fl|), for 6 = 0. 
Crosses show the result before cancellation of the 
O(a) part by the pure gauge boundary countert- 
erm, and octagons the result after the cancella- 
tion. The line shows the 0{a) part of the lattice 
artefact that is cancelled by the counterterm. 

2.1. Lattice artefacts 

The vanishing of s\ indicates the absence of 
bulk 0(a) artefacts. The non-vanishing of r[ av ^ 
reflects the presence of boundary 0(a) effects. 
These can be absorbed into a pure gauge bound- 
ary counterterm || . 

But higher order lattice artefacts, both from 
the bulk and the boundary, remain. We can study 
them for the fermionic contribution to the step 
scaling function |2|,[|] with scale factor 2. We com- 
pare the fermion contribution on the lattice (per 
continuum flavor) with its continuum limit: 



5 ltl {a/L) 



p[ a j\2L/a)-p[f(L/a) 
26o,ilog2 



(5) 



5i ) i(a/i) is shown in Figs. [I] and ^ for 9 = and 
7r/5, respectively. Its deviation from 1 at finite 
a/L is a lattice artefact. As can be seen from the 
two cases, the higher order lattice artefacts can 
depend very sensitively on the observable consid- 
ered, here the step /3-function for different values 
of the spatial boundary conditions. 
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Figure 2. Same as Fig. [I], but for 9 = ir/5. 



3. APPLICATION TO 16 FLAVOR QCD 

For n/ between 9 and 16 the perturbative two- 
loop /3-function has a second zero, that for nf = 
16 occurs at a coupling, g 2 w 0.5, where the per- 
turbative prediction might be trusted. In a recent 
study of nf — 16 flavor lattice QCD with stag- 
gered fermions we found a bulk first order tran- 
sition from a confined, chirally broken phase at 
strong coupling to a phase that appeared decon- 
fined and chirally symmetric || . Measurement of 
masses from correlation functions with pion, rho 
and nucleon quantum numbers indicated that the 
/3-function in this weak coupling phase might have 
the opposite sign than few flavor QCD, as pertur- 
bation theory predicts. However, in a deconfined 
phase it is not so clear, what is meant by "hadron 
masses" and how to relate their behavior to the 
/3-function, which is defined as the change of the 
coupling under a change of scale. 

The Schrodinger functional formalism gives a 
direct definition of a running coupling as a func- 
tion of scale, given by the size of the system. 
We have measured this coupling at several val- 
ues of the bare coupling in the "weak coupling 
phase" , and for several system sizes, in numer- 
ical simulations. These simulations have been 
done at zero quark mass, which is possible with 
the fixed boundary conditions in the Schrodinger 



Table 2 

The Schrodinger functional coupling, g 2 (L), for 
nf = 16, with 1/g 2 averaged over the choices T = 
Li 1, for various bare couplings and system sizes. 
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L = 4 


L = 6 


L = 8 


L = 12 


4.5 


5.31(15) 


3.50(16) 


3.22(14) 


2.87(12) 


4.6 


4.41(16) 


2.97(5) 


2.77(9) 




4.7 


3.92(11) 


2.78(4) 






4.8 


3.30(6) 


2.67(3) 


2.37(3) 




4.9 


3.07(7) 


2.42(3) 






5.0 


2.91(7) 


2.28(3) 







functional. This avoids having to deal with the 
anomalous dimension of the quark mass. The re- 
sults, with 1/g 2 — the observable in the sim- 
ulations (see ^) — averaged over the choices 
T = L ± 1, are listed in Table |. 

We see, from Table 0, that at fixed bare cou- 
pling, /3, the running coupling decreases with in- 
creasing system size. This is just the opposite 
behavior than that exhibited by few flavor QCD: 
the /3-function of the renormalized coupling has 
the opposite sign. On the other hand, demand- 
ing that g 2 (L/2,P') = g 2 (L,(3), i.e. that the lat- 
tice spacing at bare coupling (3' is a factor of two 
larger than at f3, we see that for n/ = 16 f3' is 
larger than /3, again, opposite the behavior of few 
flavor QCD. Hence the /3-function of the bare cou- 
pling has the opposite sign. These findings are in 
agreement with the prediction from perturbation 
theory, and confirm the results of 0. 
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